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The phase diagram of finite isospin, zero temperature QCD with the pions coupled to photons in 
a uniform external magnetic field is explored in the low field, small isospin density regime for which 
chiral perturbation theory is a valid description. For realistic pion masses, the system behaves 
as a type-II superconductor: a uniform superconducting state is formed at low-enough magnetic 
fields, a vortex state for intermediate magnetic fields and finally a normal state for large magnetic 
fields. In each these phases (including the vortex phase), the tt® remains uncondensed just as in the 
zero-external field problem. The critical magnetic field where the phase transition from the uniform 
superconducting state to a vortex state occurs was found numerically. 
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I. INTRODUCTION 

QCD in a background magnetic field is of interest due 
to its potential relevance to heavy ion collisions, neutron 
stars and perhaps even the early universe. The study 
of the QCD at zero density demonstrates a rich array 
of physical phenomena in the presence of a background 
magnetic field. The system exhibits dimensional reduc¬ 
tion which leads to stronger pairing between quarks 
and therefore enhances chiral symmetry breaking [Mg. 
Studies have also investigated how an external magnetic 
field affects the deconfinement transition temperature in 
QCD; while the magnetic field enhances chiral symmetry 
breaking, larger temperatures drive the vacuum towards 
the restoration of chiral symmetry IIIZI- Furthermore, it 
has also been suggested that for large enough magnetic 
fields, the QCD vacuum may exhibit superconductivity 
due to the condensation of p mesons. |8] 

While QCD at zero density in a magnetic field can be 
studied by both analytical and lattice methods, a more 
relevant question in the study of neutron stars relates 
to how magnetic fields affect finite density matter. This 
problem is largely unsolved except for a limited regime 
of asymptotically large baryon densities where a phase of 
color superconducting is expected to form |g. In partic¬ 
ular, the effect of magnetic fields on baryons may signifi¬ 
cantly affect the properties of magnetars—neutron stars 
with large magnetic fields nnj. However, the study of fi¬ 
nite baryon densities problems using lattice methods are 
hindered by the fermion sign problem mHig. It should 
be noted that neutron stars, possess not only a finite 
baryon density but also a finite isospin density, which 
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arises due to isospin asymmetry. The study of QCD at fi¬ 
nite isospin density (and zero baryon density) is, however, 
unencumbered by the fermion sign problem miiiiiini- 
[M| . While this system, differs from the physical system 
of relevance in neutron stars, it does give a distinct han¬ 
dle to get insights into how QCD behaves in extreme 
environments. 

Recently, finite isospin QCD in a background mag¬ 
netic field and finite temperatures was studied on the 
lattice m- The study shows the existence of the fermion 
sign problem at finite isospin QCD with a magnetic field 
even though the sign problem is absent in finite isospin 
QCD when the magnetic field is absent. This occurs due 
to the breaking of flavor symmetry of the finite isospin 
system by the external magnetic field. Note that the elec¬ 
tromagnetic charge of the up and down quarks are oppo¬ 
site in sign but different in magnitude unlike their chem¬ 
ical potentials, which are equal and opposite. The study 
circumvents the fermion sign problem by doing a Taylor 
expansion in the magnetic field. The results in that work 
suggest that at low temperatures, the system exhibits a 
uniform diamagnetic phase (i.e. negative magnetic sus¬ 
ceptibility) for isospin chemical potentials that are larger 
than the pion mass but not asymptotically large. 

In the present study, we study finite density isospin 
matter in a uniform magnetic field in a regime in which 
analytic calculations are legitimate. We work at lead¬ 
ing order in chiral perturbation theory (xPT) and con¬ 
sider the regime, where all physical parameters are much 
smaller than the typical hadronic scale (Ajjad ~ 

More specifically, the ratio 


where x can be the isospin chemical potential {pi), pion 
mass (m^n-), pion momenta (p) or V eH, where e is the 
charge of a pion and H is the external magnetic field. 
This ensures the validity of chiral perturbation theory. 
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Additionally, we assume that the baryon density is zero 
and only consider the effects of the isospin chemical po¬ 
tential. We show the existence not only of a uniform 
diamagnetic, superconducting state which persists at low 
magnetic fields but also show the existence of a topolog¬ 
ical, vortex phase. 

A system of finite isospin density was first considered 
by Son and Stephanov in Refs. [HIH], where it was seen 
that for isospin chemical potential (/x/) larger than the 
pion mass it is energetically favorable for pions to 

condense. In that work, both weak and electromagnetic 
interactions were turned off. Here we continue to ignore 
pion decay via weak interactions; this is innocuous since 
the contributions of weak interactions are small. The 
electromagnetic interactions among the pions, however 
are potentially problematic: if a system with isospin den¬ 
sity ni is confined to a region of spatial extent, R, then 
energy due to pion-pion electrostatic interactions scales 
as 


^electrostatic ^ ^ 71 jR 


( 2 ) 


where nj is the isospin density which at these low densi¬ 
ties is effectively the density of pions. The energy asso¬ 
ciated with strong interactions and the interaction with 
the external magnetic field in the regime of validity of 
xPT scales as 


^xPT 


4 /m| nj- eH\ 

Pj 




( 3 ) 


where g is a calculable dimensionless function of order 
unity. Clearly in the thermodynamic limit of i? —>■ oo, 
A-eiectrostatic dominates. Nevertheless, we will neglect it 
here even though we work in the thermodynamic limit. 

It is important to justify the neglect of the electro¬ 
static interactions between pions. Here, the calculations 
are not being used to directly describe a realistic phys¬ 
ical situation but rather to gain some insight into how 
QCD behaves by working in a particular more tractable 
regime. The issue is somewhat analogous to studies of 
infinite nuclear matter |16) . Clearly infinite nuclear mat¬ 
ter has the same difficulty—the electrostatic energy per 
particle diverges. Despite this, knowledge of the proper¬ 
ties of infinite nuclear matter neglecting the electrostatic 
energy give important insights into nuclear physics even 
in regimes where the electrostatic energy cannot be ne¬ 
glected. 

For real world parameters, there are two regimes for 
which the approximation of simultaneously neglecting 
electromagnetic interactions between pions and working 
in the thermodynamic limit appears to be appropriate. 
One is the regime in which the system is confined to a 
finite region satisfying 


1 
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J TT 
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The first condition is what is needed to justify the ther¬ 
modynamic limit when other parameters are in the typ¬ 
ical regime of validity of yPT; the second condition is 


what ensures that the electrostatic energies are small. 
Note that in the Heaviside convention used here, 1/e^ ^ 
137 so that there exists an R? which is an order of mag¬ 
nitude bigger than 1//^ while being an order of magni¬ 
tude smaller than . An alternative regime is a ther¬ 
modynamically large one in which the positive electric 
charge of the isospin matter is neutralized by a back¬ 
ground of electrons. This situation is complicated by the 
fact that these “background” electrons also couple to the 
external magnetic field and are known to exhibit super¬ 
conducting properties which might affect the behavior of 
the pions. However, the fact that electrons are so much 
lighter than the pions renders the electrons innocuous 
over the region of interest. As shown in Refs. 123) . 
the zero-angular momentum, positive parity (0^) pair¬ 
ing of relativistic electrons can exhibit either type-I or 
type-H superconductivity. Since electron masses are neg¬ 
ligible compared to pion masses, electrons in our system 
move ultra-relativistically, where the behavior of the su¬ 
perconductor is strongly type-I. The magnetic field (H®) 
at which the phase transition from a superconducting 
state to a normal state occurs is of ^(/cf^bT'c) at zero 
temperature. Here fcp is the Fermi momenta, /cb is the 
Boltzmann constant (an in our units is just 1) and Tc is 
the critical temperature below which superconductivity 
occurs. As will be shown later, the magnetic fields re¬ 
quired for which pions undergo phase changes are much 
larger and of 0{p), where /,r is the pion decay con¬ 
stant. Therefore, except for a narrow regime below 
the superconducting properties of the ultra-relativistic 
electrons do not interfere with that the pions. 


II. xPT LAGRANGIAN AT FINITE ISOSPIN 


We begin with a brief review of the finite isospin La- 
grangian at leading order in xPT that was first consid¬ 
ered by Son and Stephanov in Refs. The relevant 

regime where the Lagrangian is valid is when 1, 

where x could be the pion momenta, pion mass (jnP, 
the isospin chemical potential fj,i or y/eH, where H is 
the external magnetic field and e the pion charge. The 
effective Lagrangian has the form: 

/leff = fl Tr(D^EtD^E) + + E^) , (5) 

where E represents SU(2) matrices, ttItt the pion mass 
and the pion decay constant. The covariant deriva¬ 
tives are defined as follows: 


D^E = i9^E - i [(5^oW, S] (6) 

with the isospin chemical potential entering the La¬ 
grangian as the zeroth component. The values of the 
pion mass and pion decay constant are approximately 
135 MeV and 93 MeV respectively. In order to proceed. 
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we choose the following ansatz for the SU(2) matrix S: 


E = (crl + il^xTl + myT2 + i'KzT^) 

J-K 

= cos '(/'(cos 01 + i sin 6 *(cos ari + sin 0 x 2 )) + i sin 

m 


Note that the neutral pion is represented by and the 
positively and negatively charged pions by ± iiTy re¬ 
spectively. Plugging this ansatz into Eq. , the effective 
Lagrangian becomes: 




2 L 


cos"^ ^{sin 9{Vay -I- (V0)"^} -|- (V^/>) -I- {cos 9 cosip — 1) -f 


sin^ 9 cos^ Ip 


( 8 ) 


r 


Note that the lagrangian has been normalized such that 
CeS = 0 when 0 = 0, which is the normal QCD vacuum 
state at zero isospin or as shown in Refs. [TH [T5] . the 
vacuum at finite isospin for isospin chemical potentials 
less than or equal to the pion mass. 

From the Lagrangian density, it is straightforward to 
deduce the ground state of the system that was worked 
out in Refs. [H IH]- In order to so, we assume that 
the kinetic energy is zero and maximize the remaining 
Lagrangian. Since each of the potential energy contribu¬ 
tions are of the same sign and is a function of cosip, 
the “energy” (more specihcally % + ^ n, where H is 
the Hamiltonian, fj, is the chemical potential and n is 
the number density) is minimized for ip = 0. Then 
the Lagrangian has to be maximized with respect to 9, 
which leads to two possible solutions: either sin 0 = 0 or 
cos 9 = ^. If sin 0 = 0, then £eff = 0 but if cos0 = 

fif ’ ’ 

r2 f 2 _ 2 >,2 

then £eff = . Therefore, when the isospin 

chemical potential is greater than the pion mass, the 
condensed phase (0 7 ^ 0 ) is energetically more favorable 
compared to the normal phase (0 = 0 ). 


III. GIBBS FREE ENERGY 

Here we want to study the effect of coupling the pions 
to dynamical photons and a uniform, external magnetic 
field. The effective Lagrangian for this particular sys¬ 
tem is easily obtained by introducing photon helds and 
changing the covariant derivation of Eq. to include 
electromagnetic gauge helds. The new covariant deriva¬ 
tive is as follows: 

DyT. = dyT. - i [Sy^oni, E] - ieA^[Q, E] , (9) 

with Q being the charge matrix for the quarks and is 
dehned as: 

Q = + 2’'’3 j (10) 

where 1 is a 2 x 2 identity matrix and T 3 is the third 
Pauli matrix. The resulting effective Lagrangian using 
the dehnition E from Eq. 0 is as follows: 


Eeff = — 


1 ■ f 

4 2 


2 r _\2^ 

cos^ i/;{sin^ 0 f Va-I-eAJ + {V9^} + {VipY 


( 11 ) 


where Fij = diAj — djAi is the electromagnetic tensor. 
We have assumed here that the zeroth component of the 
four-potential A^ vanishes. 

In order to consider the thermodynamics of the h- 
nite isospin system (with photons) coupled to a uniform 
background magnetic held, it is standard to consider the 
Gibbs free energy density, which is dehned as [^ : 

g = nes-M-H ■ ( 12 ) 

Here the magnetization M is dehned as M = B — H, 
with B = W X A. As it is standard, we will assume 
that the external magnetic held, which we will label H 
only has a z-component. Since we are considering a time- 


independent system, the Hamiltonian density is given by 
the relation HeS = —^eS, where the CeS is dehned in 
Eq. (fnl). 


IV. SUPERGONDUCTIVITY 

At hnite isospin, the condensed phase of pions is a su- 
perhuid with one of the charged pions forming a massless 
mode while the neutral pion and the other charged pion is 
massive mm]. It is natural to expect then that this su- 
perhuid phase exhibits superconducting behavior in the 
presence of an external magnetic held: the system sets up 
currents that produce opposing magnetic helds to cancel 
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out the external magnetic field. For the finite isospin sys¬ 
tem in an external magnetic field, a natural question that 
arises then is the nature of superconductivity that the 
system exhibits, either type-I or type-II. The paradigm 
for studying superconductivity is Ginzburg-Landau the¬ 
ory or the Abelian Higgs model. However, the effective 
Lagrangian of Eq. (11) for the finite isospin system is 

or con- 


neither isomorphic to Ginzburg-Landau 
sequently the Abelian Higgs model. Also note that the 
finite isospin system consists of an additional degree of 
freedom (namely the neutral pion) in addition to the 
charged pions. Therefore, relative to Ginzburg-Landau, 
there is an extra degree of freedom in our problem. Fur¬ 
thermore, the effective Lagrangian for our system seems 
highly non-linear due to the presence of infinite number 
of non-linear terms present through the sine and cosine 
functions that characterize the pion fields through E de¬ 
fined in Eq. 0 - 

While the Gibbs free energy of our problem seems quite 
different from that of either Ginzburg-Landau theory or 
the Abelian Higgs model, it is still useful to borrow in¬ 
sights from the well known problems. It is well-known in 
these theories that the system exhibits type-I supercon¬ 
ductivity if: 


and type-H superconductivity if: 

^<V2X . 


(13) 


(14) 


Here ^ is the coherent length of the macroscopic state, 
and A is the penetration depth of the external magnetic 
field. The superconductivity type is determined by con¬ 
sidering the surface energy of the interface consisting of 
a half-infinite normal state and a half-infinite supercon¬ 
ducting state at the critical magnetic field, where the 
normal and superconducting states have the same Gibbs 
free energy. If the surface energy is positive at the crit¬ 
ical magnetic field, then the system makes a first order 
transition from the superconducting state to a normal. 
However, if the surface energy is negative, then an in¬ 
termediate state of Abrikosov vortices are formed such 
that the transition to the normal state from the super¬ 
conducting state becomes second order. 

While the exact mathematical results of when super¬ 
conductivity is type-I or type-H does not directly apply 
for the case of finite isospin QGD, it is still useful to com¬ 
pare the relative sizes of the coherence length and pene¬ 
tration depth in our problem. Since the coherent length 
and penetration depth are related to the masses of the 
pions and the photons, we proceed by expanding the ef¬ 
fective Lagrangian of Eq. (11) about the ground state of 

the condensed, superfluid phase 0 gs = arccos and 

^pgs = 0. In doing so, we find that the photon mass toa 
is: 


• 2 /I /• / -i ^TT 

rriA = e/^ sm O^s = QJtt [ ^ - 4 


(15) 


and the mass of the charged pion is given by 


me = y(m^AX3S^gs^-Aifcos^0g^ = /ii^ ^1 - ^ 

' (16) 

Note that here e is the pion charge, /^r is the pion decay 
constant and m^n- is the pion mass in the normal phase. 
We do not consider the mass of the neutral pion. As will 
be seen later, this is justified since the neutral pions do 
not condense. 

Since the coherent length is inversely proportional to 
the pion mass, i.e. ^ ~ ^ and the penetration depth is 
inversely proportional to the photon mass, i.e. A ^ 
it is useful to compare the masses of the pion and tfie 
photon to determine whether our system behaves as a 
type-I or type-H superconductor. 

In Fig. [l] we have a plot showing a region where the 
mass of the photon is larger than the mass of the pion (in 
orange) and a region where the photon mass is smaller 
than the pion mass (in blue). The green region is where 
the normal phase exists and is therefore, irrelevant. 
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FIG. 1: The region where itia > rng is shown in blue, the 
region where me > tua is shown in orange and the green 
region represents the normal phase i.e. pi < m,r. 

While the regions do not exactly delineate regions 
where the superconductivity is type-I or type-H, they 
do still prove some qualitative insight into the nature 
of superconductivity we might expect. For instance, it 
is reasonable to assert that in regime of the parameter 
space where ttia mg should be a type-I superconduc¬ 
tor and the regime where ttia ^ me should be a type-H 
superconductor. Therefore, there is a narrow region at 
low isospin chemical potentials and small (but unrealis¬ 
tic) pion masses (near the origin) where the system will 
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behave as a type-I (suggesting the existence of a critical 
point for small pion masses since the system goes from 
type-I to type-II with increasing isospin chemical poten¬ 
tials). However, for most of the parameter space, and 
in particular for realistic pion masses, the system will 
behave as a type-II superconductor. 


A. Type-I 


Here, we briefly consider type-I superconductivity and 
determine the critical magnetic field at which the tran¬ 
sition from the normal to the superconducting state oc¬ 
curs. First, we consider the Gibbs free energy density of 
the normal state, which occurs when 0 = 0. In this case, 
the external field completely penetrates the vacuum state 
such that the magnetic field (B) of the normal state is 
equal to the external magnetic field (ff). The Gibbs free 
energy of the system in its normal phase is 

g^ = nn-Mn-H=^H^. (17) 


B. Type-II 


Next, we consider type-II superconductivity, which as 
discussed previously is expected to occur for realistic pion 
masses. Here, we will determine the first critical mag¬ 
netic field Hci at which the transition from the uniform 
superconducting state to the vortex state occurs. We 
will also give a theoretical estimate of the second criti¬ 
cal magnetic field Hc 2 where the transition to the nor¬ 
mal state occurs. However, we do not calculate explicitly 
the surface energy, which determines whether the system 
actually behaves as a type-II superconductor. Alterna¬ 
tively, we can just compare the first critical magnetic 
field {Hci) with the critical magnetic field of Eq. (191 to 


determine that the system actually is type-II. If the first 
critical field iJci is smaller than the critical field of 
the type-I, then the Gibbs free energy of the vortex state 
is lower than the normal state and therefore, the system 
is type-II. 

The intermediate vortex state is a cylindrically sym¬ 
metric state. In order to determine the first critical 
point, we consider a single vortex state with the following 
ansatz: 


Here "Hn represents effective Hamiltonian in the normal 
phases and is the magnetization in the normal phase. 
The second equality is obtained by using Eqs. (12) and 
0 and the fact that the phase is spatially homogeneous 
with 6 = 0, '0 = 0. 

However, in the condensed phase with 6^0, the mag¬ 
netic field cannot penetrate the state. The photon fields 
readjust themselves such that the field entering the sys¬ 
tem is zero. i.e. B = 0. Therefore, the Gibbs free energy 
of the condensed phase, which is superconducting is: 


2 \2 


g, = H,-M,-H = - 


fliK-ml) 

' 2 ^ 




(18) 


where T-Lg is the Hamiltonian and Mg is the magnetiza¬ 
tion of superconducting phase. In the superconducting 
ground state 0 = arccos and 0 = 0. 

The critical magnetic field at which the phase transi¬ 
tion from the normal state to the superconducting state 
occurs when the Gibbs free energies of the normal state 
equals that of the superconducting state, i.e. = Gs- 
This gives a critical field Be, which has the following 
magnitude 


m 


- ml) 
Ml 


(19) 


a = —nep , 


( 20 ) 


where 0 is the polar angle in cylindrical coordinates, 
which goes from 0 to 27r and n is an integer. It de¬ 
termines the amount of flux that is quantized within the 
vortex. This flux ($) is given by the relation: 

$= / A-dl= f B-dA=—, (21) 

"'Coo "'Coo ® 


where the path integral is performed on the boundary 
of some region Coo, enclosing the vortex. The second 
equality arises through the use of the definition H = V x 
A and Stoke’s theorem. Note that the n that appeared 
in Eq. (20) also appears here. 


We can now write down the Gibbs free energy of a 
single vortex state with the lowest amount of flux going 
through it, i.e. n = 1. We do so using the following 
choice of gauge fields: 


A={A„ A^, A,) = (0, A^{r), 0) 


( 22 ) 


Again, we have defined the gauge choice in cylindrical 
coordinates owing to the cylindrical symmetry of the vor¬ 
tex state. The physical implication of the vortex ansatz 
is easy to understand from the definition of the electro¬ 
magnetic current, which is easily determined through the 
effective Lagrangian of Eq. (11): 


It is not surprising that the size of the critical field in¬ 
creases with increasing isospin chemical potentials. The 
density of pions that condense increases with increasing 
chemical potentials, which in turn means that larger cur¬ 
rents are generated and the superconducting state per¬ 
sists at larger external magnetic fields. 


j = ^^ = qcos^ 0sin0 ^Va -I- eAj . (23) 

Noting that the flux quantization condition implies that: 

lim A^{r) = — , (24) 

r—¥oo CT 
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which along with the definition of the electromagnetic 
current in Eq. (23), implies that the current at the bound¬ 
aries of the vortex vanishes. 


Using the flux quantization condition then, the Gibbs 
free energy of a single vortex state with n = 1 is: 


^vortex — Qt Qs 


= X - 


1 f 1 d{rA^) \ /■ 


2 V r dr 


cos^ '(/'{sin^ 


„2 f2 

— rri^f^ cos 9 cos ^-sin^ 9 cos^ 


Mt = (B, - H) 

'ld{rA^) 


Bt = z 


r dr 


where is the Gibbs free energy of the uniform super¬ 
conducting state, which was calculated in Eq. (18). 

Since it is not possible to And the solutions of the vor¬ 
tex analytically by solving for the gauge held A and the 
pion held implicitly represented by 9 and a, we will pro¬ 
ceed to numerically And the solution of the vortex state. 
An important question arises in the context of solving 
for the vortex state. It pertains to the fact that in chi¬ 
ral perturbation theory, in addition to the charged pions 
there are also neutral pions. These neutral pions are un¬ 
condensed as was shown in Section II based on results of 



^-MfH 


(25) 


Refs. [TilITK] . It is important to establish that that this 
remains the case in finite magnetic fields. 

Formally, checking that the neutral pion fields indeed 
do not condense i.e. '0 = 0 , amounts to finding solutions 
for charged pions under the assumption that the neutral 
pions do not condense and checking that the solutions are 
stable against fluctuations of the pion field. This can be 
done using the equation of motion for 0 fields. Assuming 
radial symmetry for this field, the equation of motion is 
as follows: 


9t0(r) — V"^0(r) = I 2 cos 0 sin 0sim 0 I-h eA^ 


(26) 


Expanding 0 around its ground state value, i.e. 0g.s. = Re^„ exp (zE„t)50„(r), we obtain the following equa- 

0, such that 0 = 0g.s. -I- Sip and assuming Sip = tion: 


- V'^Sipn + fvSipn = ElSipn 

/v = fsin^ 9{r) + eA 0 (r)^ -f 


/ d9{r) 
dr 


^ - ml fl cos 9 {r) 


Bill sin^ S{r) , 


(27) 


where /v is obtained by plugging the vortex solutions 
found by assuming that the neutral pions do not con¬ 
dense. If none of the the eigenvalues of the above equa¬ 
tion are negative, i.e they satisfy El > 0 for all n, the vor¬ 
tex is stable against neutral pion condensation; if < 0 


for some n, then the vortex is unstable. However, it is 
difficult to explicitly determine the sign of even just the 
smallest valued El since we can only numerically solve 
for the vortex solution, which in turn affects the exact na¬ 
ture of /v. Therefore, we will proceed by minimizing the 
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free energy of the vortex state with respect to all the de¬ 
grees of freedom including the gauge field A and the pion 
fields, 6 and ip. In doing so we find that it is energetically 
favorable for the neutral pion to remain uncondensed. 


C. Numerical Results 


In this subsection, we solve for the vortex solutions of 
the system directly by minimizing the Gibbs free energy 
of Eq. (25) and determine the first critical field. For the 
purposes of this numerical work, we will use the pion 
decay constant, /^, to set the scale in the problem. We 
will introduce the following change of variables: 



r 


h 

A(f, 


Ur 

U 

Ml 

u 

Acj, 

77 


(28) 



FIG. 2: The plot shows the density of the charged pion as a 
function of radial distance from the center of the vortex, which 
is at the origin. The result is for m,r = 1.5 and /tj = 2.0. 


The numerical work is performed by discretizing the 
Hamiltonian and minimizing the energy per unit length, 
which is obtained from the discretized Hamiltonian 
through a numerical integration. The numerical error 
is dominated by the discretization of the Hamiltonian 
and the error that arises in the energy per unit length 
is of O (a^), where a is the size of each cell. The error 

arising from the numerical integration itself is O ( ) j 


where N is the number of discretized points and Af is the 
length of the integration region. Therefore, the numeri¬ 
cal estimation of the critical magnetic field of Eq. (291, 

Mo 


in particular Hd = —jU has an error approximately of 


0{a^). For our calculation the maximum value of Ar is 


FIG. 3: The plot shows the magnetic field i.e. B = V x A 
as a function of radial distance from the center of the vortex. 
The result is for m,r = 1-5 and jh = 2.0. 


22 and N was chosen to be 160. Therefore, the numerical 
errors in the results we present are relatively small. 

For realistic pions masses, m,r = 1.5 and we solve for 
single vortex solutions up to an isospin chemical poten¬ 
tial of Mi = 2.5. These values are well within the regime 
of validity of chiral perturbation theory, which was dis¬ 
cussed earlier. In Figs.[^and[^ we plot the radial profile 
of a single vortex solution. We omit the plot for ipU) 
since it remains uniformly zero. 

In generating the solutions, we assumed that e = 1. 
Note that at the center of the vortex both 0(r) and A^(r) 
must vanish to prevent a singularity from occurring at the 
origin. (This can be easily understood from the equations 
of motion written in cylindrical coordinates.) Addition¬ 
ally, far away from the vortex 0 assumes its ground state 
value for a uniform superconducting state and the mag¬ 
netic field B = V X A also vanishes. 



FIG. 4: The plot shows the critical magnetic field for the 
formation of the vortex state for realistic pion mass: rtiTr = 1.5 
as a function of isospin density. 






































We also determine numerically the first critical mag¬ 
netic field {Hci) at which the transition from a uni¬ 
form superconducting state, which exists at low magnetic 
fields, to a state with a single vortex. It is a standard re¬ 
sult iniiE3 that the critical field is given by the relation: 

i?ci = ’ (29) 

5>o 

where E^oitex/L is the energy per unit length of the vor¬ 
tex state, which is determined by integrating the Hamil¬ 
tonian of the vortex over the cross section of the vortex, 
and $0 is a single quanta of flux passing through the vor¬ 
tex. The size of this flux is ^ with e being the charge of 
the pion. The relation itself is easily determined by us¬ 
ing the fact that at the critical held Hd, the Gibbs free 
energy of the superconducting state is equal to the Gibbs 
free energy of the single vortex state, i.e. = ^vortex- 
The hrst critical magnetic held is presented in Fig. 
The size of critical held increases with increasing isospin 
chemical potential. This is expected since larger isospin 
chemical potential results in a larger density of condensed 
pions in the condensed phase, which further means that 
larger currents to oppose the external magnetic held can 
be generated at the boundaries. 


D. Second critical field 


lap such that the average spacing between the vortices 
of 0{mg^). Each vortex, however, continues to car ry a 
single quanta of hux $o = Then, using Eq. (21), an 
estimate for the second critical point is given by 


Hc2 ~ ^omi , 


(30) 


where mg is given in Eq. (16). 


V. FINAL COMMENTS 


In this paper, we have investigated hnite isospin QGD 
in an external magnetic held at lowest order in chiral per¬ 
turbation theory. We have shown here that the system 
forms a type-II superconductor for realistic pion masses 
and possibly also a type-I superconductor for small pion 
masses. It will be interesting to see if the vortex phases 
that have been observed here are also seen in lattice cal¬ 
culations. It is noteworthy that lattice simulations of h- 
nite isospin system in a magnetic held have already been 
carried out [^. We note however, that the setup of 
this calculation apparently excluded the physics associ¬ 
ated with type-II superconductivity, for which vortices 
form and thereby alter the magnetic helds, i.e. the back- 
reaction of the pion held dynamics on the B held appears 
to be absent. 


While it is possible to numerically determine the lower 
critical magnetic held, where the transition from a uni¬ 
form superconducting state to a vortex state occurs, it 
is generally harder to determine the upper critical point 
from which the transition from a vortex state to a normal 
state occurs. As the magnetic held increases past the 
hrst critical point, the density of vortices also steadily 
increase forming Abrikosov lattices. Glose to the upper 
critical point, the cores of the vortices (the size of which is 
determined by the charged pion mass mg) begin to over- 
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